Phase diagram of localization in a magnetic field 
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The phase diagram of localization is numerically calcu- 
lated for a three-dimensional disordered system in the pres- 
ence of a magnetic field using the Peierls substitution. The 
mobility-edge trajectory shifts in the energy-disorder space 
when increasing the field. In the band center, localized states 
near the phase boundary become delocalized. The obtained 
field dependence of the critical disorder is in agreement with 
a power-law behavior expected from scaling theory. Close to 
the tail of the band the magnetic field causes localization of 
extended states. 



Thaugh being intensively investigated since the late 
1950sEl, the problem of the disorder-induced metal- 
insulator transition (MIT) of non-interacting electrons 
in three dimensions (3D) can be considered as being still 
unsolved. There is a controversial discussion, whether 
the critical behavior at the Anderson transition (AT) can 
be classified with respect to the fundamental symmetry 
of the Hamiltonian as proposed within field theoretical 
descriptionscrH. One would expect that critical proper- 
ties governed by the exponent of the localization length 
u were altered when the universality class was changed. 
For instance, by applying an external magnetic field to a 
system with spin rotation invariance, time reversal sym- 
metry is broken and the universality class changes from 
orthogonal to unitary. On the basis of considerably re- 
duced error bars recent numerical calculationgj seem to 
provide evidence that v is sensitivei-ljo symmetry break- 
ing in contrast to previous resultsElH. Therefore, it is 
particularly interesting to investigate the importance of 
the magnetic field for the localization mechanism that 
drives the MIT. 

In this paper we report results on the influence of a 
homogeneous magnetic field on the phase diagram of lo- 
calization of the disordered Anderson model (AM). Pre- 
vious findings concerned the phase diagram without mag- 
netic field, including both exact analytieal solutions for 
the infinite Cayley tree (Bejt|hfi-j£uttice)l3El and numerical 
studies for a real 3D latticeOu'ElEJ. Based on numerical 
results obtained by the transfer-matrix approach some 
conclusions on the nature of the wave furictions and lo- 
calization properties have been drawntfcl: In the band 
center, electrons are localized due to quantum interfer- 
ence, while outside the unperturbed band electrons be- 
come localized below a certain value of the density of 
states (DOS)lj. Here, we consider the magnetic field as 



an addiiional continuous parameter which can induce the 
MITOIIj. By using the transfer-matrix method we cal- 
culate a complete mobility-edge trajectory in the energy- 
disorder space for a finite magnetic field. We find two 
regimes with entirely different behavior with respect to 
the magnetic-field-driven MIT. In the band center, lo- 
calized states near the zero-field phase boundary become 
extended when applying a magnetic field. In slightly dis- 
ordered metallic systems for states close to the mobility- 
edge the field effect is opposite. _ 

We use the AM with diagonal disordertl aad Peierls 
phase factors in the hopping matrix elementsEJ, 



if = ^6r|r)(r|+^tr.r+A|r)(r + A| 



(1) 



r,A 



Energies are measured in units of the modulus of the 
hopping matrix elements tr,r-i-A i lengths in units of the 
lattice constant a = 1. The states |r) are associated 
with the sites of a simple cubic lattice. The site energies 
Er are distributed uniformly at random between —W/2 
and W/2. Only hopping between nearest neighbors r 
and r + A is taken into account. The hopping matrix 
elements 



r,r+A 



1, 



A e {±ej, 
A G {±ea;,±ej. 



(2) 



describe a system with a homogeneous magnetic field B 
in the x direction and the Peierls phase a = eB/hcJi& 
the number of flux quanta $o = hc/e per unit cellllj. 
Here, the Landau gauge with the vector potential A = 
(0, —Bz,Q) is chosen. The Hamiltonian (|l|) is symmetric 
and periodic with respect to a with a period a = 1 ($0 
periodicity). For a quasi- ID, bar-shaped system of cross 
section M x M the Schrodinger equation H\^) = E]"^) 
with 1^)1-^ X]'^r|r) can be written as a transfer matrix 
equationij = T^u^-i, where the vector 



(fli 



, aMMz 



(3) 



contains the coefficients of the planes z and z + 1 and 



I 



I 



(4) 



is the 2M^ x 2M^ transfer matrix. Each matrix is 
of the order IvP x Af^ and contains the elements, of the 
Hamiltonian in the z plane as presented in Ref.tl. The 
matrices I connect successive xy planes. The matrix is 



1 



a function of the energy E, the magnetic field parameter 
a, and the random site energies on the slice z. 

In discufdered quasi- ID systems all eigenstates are 
localizedllj. The wave function of a localized state de- 
cays exponentially with a localization length. The largest 
localization length Am can be identified as 7^^, the in- 
verse of the smallest positive Lyanunov exponent (LE) 
7 = min{x(uo),x > 0} defined bj£3 



x(uo) 



lim — In 

L — *oo L/ 



■Tiuol 



(5) 



where Uq is the initial vector for 2^ 0. 

By using the Benettin algorithniEj-we extracted 7 from 
the approximate spectrum of LEsllj. The localization 
length Am is a function of the cross section M of the 
bar, the electronic energy E, the number of flux quanta 
per unit cell a, and the disorder W . 

The critical behavior of Am near the MIT can be de- 
termined rumierically by establishing the one-parameter 
scaling lawLj 



XM(E,W,a) 
M 



f 



Xo,{E,W,a) 
M 



(6) 



where Aoo = hniAf^oo Am is the localization length of the 
3D system in the thermodynamic limit. The phase dia- 
gram of localization for a given a describes the MIT in 
the {E, W) plane. The mobility-edge trajectory E^iyV, a) 
is determined by using the property that at the critical 
point the reduced localization length of a quasi- ID sys- 
tem Aaz/AjT is independent of the cross section M, 



\M{Ec,W,a) 
M 



= const. 



(7) 



As an example, Fig.|^(a) shows the inverse of the re- 
duced localization length M/Aa/ as a function of the 
disorder W , at the band center E — Q. The Peierls 
phase was chosen to be a = 0.25, at which the mag- 
netic field has the strongest effect and the system belongs 
to the unitary universality class. The data were calcu- 
lated for different M with a statistical accuracy of 0.25%. 
The MIT is indicated by the common crossing point at 
Wc{E = 0,a = 0.25) = 18.35 ± 0.11. The sign of the 
size effect on Am changes when increasing the disorder 
from the metallic {W < Wc) to the insulating regime 
{W > Wc). In the same way we locate the critical dis- 
orders for energies wp to E = 7. By expanding Eq. (|^) 
around the critical point at a given energy and fitting the 
numerical data to the linearized form 



M 



= Ae + A(W -Wc)M^/'', 



(8) 



one can extract the critical exponent ly and the disorder 
dependence of the|-.localiZ|ation length Aoo, as has been 
performed in Refs.B andtZl. We have found = 0.568 ± 
0.076 for a = and A^ = 0.5|64 ± 0.027 for a = 0.25 
consistent with data from Ref.Q. 



For energies beyond the unperturbed band \E\ > 6, 
it is convenient to determine the critical points from the 
energy dependence of the localization length Am at fixed 
W. In Fig.|l|(b), the localization length is shown as a 
function of the energy for various M. The statistical ac- 
curacy is 1%. The intersect of the lines signalizes critical 
behavior and yields the position of the mobility edge. At 
W = 12, for example, Ec{a = 0.25) = 7.22 ±0.14. Other 
critical points Ec{W,a) were determined similarly. 

Combining all the data obtained by the finite-size scal- 
ing analysis as described above, one can construct the 
entire phase diagram of the AM in the parameter space 
spanned by E and W. Figure ^ shows the mobility-edge 
trajectory with and without a magnetic field. For zero 
field the states in region I are extended whereas the states 
in II are localized. The regimes I (metallic) and II (in- 
sulating) are separated by the mobility-edge trajectory 
(solid line). The latter is modified when a magnetic field 
is applied. For a = 0.25 states in region III become 
extended and the phase boundary shifts to higher val- 
ues of critical disorder. It is known from the theory 
of weak localizationEJ that in the presence of a mag- 
netic field coherent time-reversed paths are eliminated 
and hence backscattering is suppressed. Without a mag- 
netic field, states in III are localized by quantum interfer- 
ence effects. Ift-tius case a weak magnetic field leads to a 
delocalizationEJ^EJ. As a consequence, a stronger disorder 
is required in order to localize these states (see Fig. 2). 
Thus, insulating systems with W slightly larger than the 
critical disorder at zero field [W > Wc{a — 0)] undergo a 
transition to a metal when the magnetic field is applied. 
This behavior is consistent with the raachanism of the 
field-induced MIT proposed by ShapiroEJ. 

The increase of the critical disorder Wc with the mag- 
netic field a at ii^ = is shown in the inset of Fig. |^. 
For small a the field dependence can be described by the 
following relation: 



Wcia)-Wc{0)^iEcia) 



oc a 



E=Ea{0} 



(9) 



which 



obtained previously using the scaling 



approachta. Our data for Wc{a) are in agreement with 
this power law wherej-|thc critical exponent w 1.4 is 
taken from Refs.cl andQ. 

In region IV, the trajectory moves into the metallic 
phase in contrast to region III. Here, states that are ex- 
tended for a = become localized for a = 0.25. This can- 
not be explained by the interference mechanism discussed 
above. A qualitative understanding can be achieved by 
considering the motion of the trajectory of slightly dis- 
ordered systems when applying a magnetic field. In the 
limit W ^ the mobility edge merges with the band 
edge of an ordered system. The zero-field value of the 
band edge is Ei,{a = 0) = 6. In a magnetic field, the 
band edge as a function of a can be cpJculated numer- 
ically by using Hofstadter's algorithmlla generalized to 



2 



3D. The band shrinks if a magnetic field is apphed and 
varies with a period of one flux quantum per unit ceU, 
a = 1 ($0 periodicity). For example, for a — 0.25 the 
unperturbed band edge shifts to -Ef,(a — 0.25) w 4.8 (Fig. 
H). In slightly disordered systems this behavior persists. 
The band edge of such a system moves to a lower value 
Eb{W,a) < Eb{W,0), indicating that the DOS for ener- 
gies near the zero-field mobility edge decreases dramati- 
cally with increasing a. Below a certain value of the DOS 
an MIT is inducedQ. In region IV, this causes a shift of 
the mobility edge to lower energies EJW, a) < Ec{W, 0) 
when applying a magnetic field (Fig. 0). 

In fact, the shift of the phase boundary is due to the 
combination of both types of field effect mentioned above 
which compete with each other. Close to the band cen- 
ter (III), where the DOS changes negligibly with E, the 
interference effect dominates. On the other hand, in re- 
gion IV, close to the band tails, the DOS effect is much 
stronger than the interference effect. Here, the behavior 
of Ec{a) is determined mainly by the energy and field de- 
pendence of the DOS. For example, for a = 0.25 we find 
an intersect of the two phase trajectories at W* w 14.5 
and E* « 7.8, where the two effects are of the same order 
of magnitude. We believe that a similar field dependence 
of the phase diagram is also valid for a Gaussian distri- 
bution of on-site energies Cr- 

We now again concentrate on the band center E ~ Q 
in order to investigate the field dependence of the crit- 
ical parameters. Figure |3| shows the localization length 
Xi\i of quasi-lD systems as a function of magnetic field 
for various W corresponding to the delocalized, critical, 
and localized regime in zero field, respectively. The data 
were calculated with a precision of 0.25%. For a differ- 
ent from half integers, \m{E = 0) is larger than without 
a magnetic field. This field-induced enhancement of the 
localization length is directly related to the shift of the 
phase boundary in region III. One sees in Fig. 3 that for 
E = 0,W > 16.5, and M = 4 the localization lengths 
Xm coincide for both a and a — 0.5 within the sta- 
tistical uncertainties. They vary with a period of half a 
flux quantum per unit cell ($0/2 periodicity). This was 
also checked for other widths of A/ = 5, ... , 12 by com- 
paring Xm{E = 0, a = 0) with Xm{E = 0, a = 0.5). 
In the insulating region our results are consistent with a 
periodically varying transition amplitude of strongly lo- 
calized electrons as derived analytically with a directed 
path method by Lin and NoriE3. However, we find that 
the $0/2 periodicity does not persist for i? ^ as shown 
in Fig. 3. Thus the former seems to be an intrinsic prop- 
erty of the band center, around wkich the on-site energies 
tr are distributed symmetricallyEil. Furthermore, for the 
investigated systems of width M = 4, . . . , 12 we observe 
that AA/(a = 0.5) < XM{a — 0), though for both values 
of a the Hamiltonian belongs to the orthogonal univer- 
sality class. Assuming that Ac does not change in this 
case, one obtains that for a — 0.5 and E ^ {) the MIT 
should unexpectedly occur at a lower disorder than in 
the zero field case, Wc{a — 0.5) < Wc{a ~ 0). 



The maximum-entropy ansatz predicts the increase of 
the localization length in quasi- ID systems by a universal 
factor of 2, when breaking the time reversal symmetrj£3. 
Since in the present calculations the parameters are be- 
yond the range of validity of this universal relation, it 
would be desirable to extend the investigations to this 
regime for checking the prediction. Another interesting 
problem is to stud#,the disappearance of this relation, as 
was argued in RefBj, for d > 2, when extrapolating to 3D 
systems by scaling the cross section M of the quasi- ID 
bar. 

In conclusion, we have numerically calculated the 
phase diagram of localization in 3D disordered systems 
in the presence of a magnetic field. Comparing the ob- 
tained diagram with the zero-field result, we identify two 
regimes with different magnetic field dependence of the 
phase boundary. In the band center, the phase bound- 
ary is shifted towards higher values of critical disorder, 
so that the metallic phase is broadened. This is mainly 
due to the suppression of the interference of time-reversed 
paths by a magnetic field, leading to the delocalization 
of electron states. On the other hand, close to the band 
tails the location of the MIT in slightly disordered sys- 
tems (W^ ^ 0) is dominated by the field dependence of 
the DOS. Here, the mobility edge shifts towards smaller 
energies, thus diminishing the metallic phase. Our nu- 
merical findings for small fields show that the behavior 
of the critical [disorder is consistent with predictions by 
scaling theoryts. Thus, we have shown that the phase 
diagram of localization is infiuenced by an external per- 
turbation which breaks the time reversal invariance. 
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FIG. 1. Disorder and energy dependence of the reduced 
localization length Xm/M of quasi- ID disordered systems 
for various cross sections M in a magnetic field a — 0.25. 

a) M/Xm vs. disorder W at the band center E — 0. 

b) \n{XM /M) vs. energy E at W = 12 outside the unper- 
turbed band. Lines are polynomial interpolations. 

FIG. 2. The phase diagram of the Anderson model of lo- 
calization. Full dots show the critical points {Ec, Wc} for a 
magnetic field a — 0.25. The error bars indicate the statis- 
tical uncertainties due to the finite length of the quasi-lD 
systems. The OBen circles show numerical results for q = 
taken from Ref.B. The full and the dotted line represent the 
mobility edge trajectory for a — and a — 0.25 respectively. 
Inset: critical disorder Wc vs. magnetic field a in the band 
center E = 0. The-curve shows the power-law behavior from 
scaling argumentslij. 

FIG. 3. Localization length Am of a quasi-lD disordered 
system with a cross section M = 4 for various disorders W 
corresponding to different regimes (extended states, critical 
region and localized states) as a function of the number of 
flux quanta per unit cell a for E — (full dots) and E = 4 
(open circles). The continuous lines are fits by <l>o/2-periodic 
functions. 
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fig. 1, "Phase diagram of localization in a magnetic field", Drose et al. 




fig. 2, "Phase diagram of localization in a magnetic field", Drose et al. 




fig. 3, "Phase diagram of localization in a magnetic field", Drose et al. 



